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ABSTRACT 



Water-type solar heating systems have received much 
attention in the literature and several good computer simulation 
models are available. Air-type space heating units have 
advantages which may make them more suited to a particular 
need than the water- type system. 

This computer simulation offers a model of an air-type 
solar heating system with which one can determine the monthly 
system output. All significant system parameters are easily 
varied as is the storage bed size and material. The simulation 
accepts all air-type collector designs. The residence is 
simply modeled and estimating its heat load is reduced to 
calculating the product of floor area times a constant. 
Insolation data input is based on monthly averaging- techniques . 

The quantities varied in the present application of the 
model were collector size and storage volume. 
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NOMENCLATURE 


Numbers 


in (. ) refer to the equation (s) in which 
the symbol appears. 


A (12,14,15) 


cross-sectional area of pebble bed storage 
unit m 2 


A c 


overall collector area m 2 


A p (13) 


total surface area of pebble bed storage 
unit m 2 


A xs (9) 


area of collector normal to air flow m 2 


b 


subscript refering to the storage bed 
material 


c b (13,15) 


specific heat of pebbles in storage bed; 
assumed constant at 0.837 KJ/Kg - °C 


c f (11-13,15) 


specific heat of air; assumed constant 
at 1.012 KJ/Kg - °C 


D (14) 


diameter of pebbles, all assumed equal m 


D h (7,9) 


characteristic length of collector 
(hydraulic diameter) for use in determining 
Reynold's number; equal to twice the 
spacing between the plates m 


d (23) 


dust factor 


F' (11) 


collector efficiency factor [1] 


F r (10,11) 


collector heat removal factor [1] 


f 


subscript refering to the fluid (air) 


G (14) 


flow rate per unit cross-sectional area 
of storage bed Kg/s - m 2 


G c (ID 


flow rate per unit of collector area 
Kg/s - m 2 


H (16,23) 


monthly average daily total radiation on a 
horizontal surface KJ/m 2 - day 
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h 1 (7) 


convective heat transfer coefficient for 
collector absorber plate W/m 2 - °C 


h 2 (7) 


convective heat transfer coefficient for 
collector back plate W/m 2 - °C 


h r (6) 


radiation heat transfer coefficient between 
the two air duct surfaces of collector 
(sides neglected) W/m 2 - °C 


h v (12,14,15) 


volumetric heat transfer coefficient for 
storage bed material W/m 3 - °C 


K d (16,23) 


ratio of the monthly average daily diffuse 
radiation to the daily extra-terrestrial 
radiation 


K t (16,23) 


ratio of the monthly average daily total 
radiation to the daily extraterrestrial 
radiation 


k air (7) 


thermal conductivity of air, taken to be 
constant at 0.029 W/m - °C 


k ins 


thermal conductivity of insolation used 
under collector, taken to be constant at 
0.043 W/m - °C 


L (12,15) 


length of storage bed m 


m (12,13) 


subscript referring to the mth segment 
of the storage bed 


m (9,12,13,15) 


mass flow rate of air Kg/s 


N (1,25) : (20) 


number of glass covers: day of year measured 
from January 1 


Nu (7,8) 


Nusselt Number 


n (12,13,15) 


number of segments into which the storage 
bed is considered to be divided for 
application of numerical technique 


n g (27) 


refractive index of glass 
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Q u (10) 


hourly rate of useful thermal energy 
collection per unit area of collector 
surface W/m 2 


Rjj (16,17,23) 


ratio of radiation on a tilted surface 
to that on a horizontal surface 


R p (8,9) 

D,H 


Reynold's Number for collector air flow 
based on hydraulic diameter 


r d (16,21,23) 


ratio of hourly to daily diffuse radiation 


r t 


ratio of hourly to daily total radiation 


S (19,23) 


hourly flux of solar radiation absorbed 
by collection surface W/m 2 


SI (16) 


hourly flux of solar radiation striking 
a tilted surface W/m 2 


S (23) 


shading factor 


T amb (1 ' 10) 


ambient temperance °C 


T c,i (10) 


inlet air temperature to the collector °C 


T b (12,13) 
m 


temperature of pebbles contained in storage 
segment m 


T f (12,13) 
m 


temperature of air entering storage segment 
m °C 


T f (12,13) 

r m+ll 


temperature of air leaving storage segment 
m+1 °C 


T f ,m 


mean fluid temperature in collector 


Tp (1) 


mean absorber plate temperature during 
operation of collector °K 


T s (13) 


temperature of the surroundings containing 
the storage bed °C 


T (6) 


mean temperature for radiation between 
the absorber plate and back plate of 
collector, assumed constant at 340 °K 
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T x (6} 


mean temperature of absorber plate °K 


T 2 (6) 


mean temperature of back plate °K 


^ns ^ 


thickness of insulation behind back plate m 


U b (5) 


overall collector back loss coefficient 
W/m 2 - °C 


U L (10,11) 


overall collector loss coefficient 
W/m 2 - °C 


U LS (13 > 


overall storage loss coefficient W/m 2 - °C 


U t (45°) (1,4) 


overall collector top loss coefficient 
for a unit tilted 45 from the horizontal 
W/m 2 - °C 


U t (B) (4) 


overall collector top loss coefficient 
for a unit tilted degrees from the 

horizontal W/m 2 - °C 


V (3) 


wind velocity, assumed constant at 5 m/s 


a 


solar altitude 


a(0 t ) (24,28) 


directional absorptivity of collector plate 


8 (4,16,23) 


collector tilt from horizontal 


Y (18) 


surface azimuthal angle (due soth being 
zero and east positive) 


6 (18-20,22) 


declination 


e (13,15) 


void ratio of storage bed 


e g 


emissivity of glass cover plates 


£ P U) 


emissivity of absorber plate 


e 1 (6) 


emissivity of absorber plate 


e 2 (6) 


emissivity of back plate 


0 t (17,18,24-27) 


angle of incidence of direct radiation 
measured from the normal to the collector 
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0 2 (26,27) 


angle of refraction 


y (9) 


absolute viscosity Kg/m-s 


p (13,15) 


density of pebble bed material Kg/m 3 


P gr (16,23) 
p(9 t ) (25,26) 


diffuse ground reflectivity 

reflectivity of a single air-glass inter 
face for direct radiation as a function 
of <f> t 


p (60 ) (24) 


reflectivity of a single air-glass 
interface for diffuse radiation 


T (12,13,15) 


time s 




effective transmittance - absorptivity 
product for direct (beam) radiation 
and diffuse radiation, respectively 


T(6 t ) (24) 


directional transmittance of glass as a 
function of 0 t 


<f> (18,19,22) 


latitude (north positive) 


a) (18,19,21) 


hour angle (solar noon zero and morning 
positive) 


u) s (21,22) 


hour angle corresponding to sunrise 
(and sunset since the day is assumed 
symmetric about solar noon) 


a = 


= 5. 6693 (10) ” 8 W/m 2 - °K 4 
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I . INTRODUCTION 



This study was undertaken for the purpose of developing 
a computer simulation of an air-type residential solar 
space-heating model from which to determine the optimum 
collector and storage size for a given heat load. Furthermore, 
the need for comparison between this and liquid-type systems 
dictated that the model developed be adaptable to changes 
in various significant parameters and components. 

To date the majority of published computer simulation 
models, as well as actual installations, has been of the 
" liquid-type" , that is, a solar heating system wherein 
water or a water-antifreeze solution is the working fluid [7]. 
The most noteable exception to this trend has been the 
simulations and installations of Dr. G. Lof. 

Liquid-type systems are favored [7] for commercial 
installations where a large amount of energy must be 
distributed. In such a case the storage tank, piping and 
insulation required occupies far less valuable space than 
would the storage and ductwork necessary to transport the 
same amount of heat using air at reasonable rates of flow. 

Also, the basic components (heat exchangers, pumps and 
controls) associated with the liquid-type system become 
proportionally less expensive as the size of the system 
they serve grows. 
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On the scale of a single- family dwelling however, it 
appears a trade-off becomes practicable between liquid and 
air systems depending upon the geographic location, nature 
of an existing conventional system (for retrofit solar 
installations) and extent of solar energy utilization desired, 
among other considerations. Some salient trade-off 
considerations of the two types of systems currently most 
widely employed are summarized in Table 1. Though not 
exhaustive, this indicates generally some of the more 
significant differences between the systems. 

These points serve only as a first-cut at making an 
overall system evaluation because various design strategies 
may be employed to lessen the disadvantages of the liquid 
system or extend the utilization of the air system. 

For example, if in a water-type system the space 
heat is supplied by blowing air over a heat exchanger 
directly into the room, instead of relying on natural 
convection from baseboards or radiators, lower water 
temperatures may be acceptable. Solar system efficiency 
increases significantly when the thermal energy is absorbed 
at the lowest useful temperature. This heat transfer 
can be accomplished at the storage container (e.g., in a 
rock bed or finned annulus surrounding the tank) with 
standard ducting distributing the heated air throughout the 
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TABLE 1: Comparison of relative merits of linuid and air-type 

solar heating systems. “Liquid*' refers water or 
water-antifreeze solution as the thermal energy 
transport medium coupled to a water storage tank. 
"Air" refers to air as the transport medium with a 
rock bed for storage. 

The base for comparison is space heating with 
additional abilities treated as advantages over this. 



LIQUID 

ADVANTAGES 

1. smaller storage volume per 
unit of heat stored 

2. transport network occupies 
smaller space 

3. directly adaptable for use 
with absorption A/C 
equipment 

4. directly adaptable to 
domestic hot water supply 

5. adaptable to either forced- 
hot-water or forced-hot-air 

6. if chilled- water air-condi- 
tioning is supplied, one 
distribution system serves 
for both heating and cooling 



DISADVANTAGES 

1. more complex controls, heat 
exchangers required in most 
applications, maintainance 
more costly 

2. Less efficient energy collection 
at temperatures required for 
forced-hot-water baseboard 
heating (140-190OP) 

3. may require some treatment of 
transport liquid, such as: 
antifreeze, antiscale, 
anticorrosive or biocide 

4. possibility of expensive damage 
in dwelling due to collector 
leakage, or to collector duo to 
freeze-up 

5. water may bo scarce commodity; 
antifreeze (oil by-product) has 
uncertain future cost and 
availability 

6. where central air conditioning 
is supplied, forced-hot-water 

heat requires parallel distribution 
system (pipes) for space heating 



AIR 

ADVANTAGES 

1. simple operation, longer lifetime, 
loss maintainance, less capital 
and installation cost 

2. since the transport fluid is 
identical to the ultimate delivery 
medium, no heat exchanger, per sa. 
Is necessary 

3. system can be adapted, though not 
readily, to cooling of dwelling* 

4. where central air conditioning is 
supplied, one distribution system 
serves for both heating and cooling 
(ductwork sized for the rua®s f 
rates associated with A/C lends 
itself to collection of heat at the 
lowest temperature) 

5. air is free, rocks are plentiful 
and leaks are neither dangerous 
nor damaging 



DISADVANTAGES 

1. not readily adaptable (possibly not 
feasible) to providing source of 
heat for conventional air 
conditioning equipment 

2. where (finned tube) heat exchangers 
are used, they must be larger than 
that required in a liquid-system 

to transfer the same quantity of: 
energy 

3. ducting and storage reauire more 
space for the same energy transport 
or storage, respectively 



*e.g., cool storage at night for use 
during the day. 
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dwelling. Alternatively the heated water can be piped to 
those locations where space heat is required and air blown 
over a coil-type heat exchanger into the space. The latter 
scheme might be advantageous where space is at a premium 
or where the distribution network is extensive. 

On the other hand, by placing a water coil within the 
rock bed of the air-type system or girdling the bed with 
such a coil, preheated or fully heated domestic hot water 
can be provided. Alternatively, by cooling the rock bed 
by night with ambient air, the dwelling can be cooled 
during the day using the storage as a heat sink. The materials 
of the rock bed being rather low in cost, it might be 
worthwhile to provide a separate , uninsulated storage bed 
a meter or so underground to be cooled at night in the 
above manner while also losing heat to the ground; this 
unit would be exclusively for summer air conditioning 
and sized accordingly. The main storage unit, well insulated 
and connected to the collector, would be used to heat 
domestic hot water (as described above) and would be 
available, as well, to provide space heating as needed 
on cool summer nights or days. This latter bed would be 
sized for winter space heating. 

As can be seen, hybridization and innovation can be 



used to great advantage, especially in adapting to an 
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already existing conventional system, but makes a clear-cut 
determination of the "best" system for a given situation 
very difficult. Regional weather characteristics and labor 
rates, the cost in dollars of one system versus another 
(including various levels of component refinement in each 
case) and the availability of materials (e.g. a reliable 
source of water or antifreeze) only complicate further a 
fundamentally complicated problem in thermodynamics . To 
this end, "canned" computer simulations appear to offer the 
most direct and cost-effective manner of gleaning from 
the myriad of interrelationships the optimum design for a 
given dwelling and location. 
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II. DESCRIPTION OF THE SOLAR HEATING 
SYSTEM SIMULATED 



A schematic diagram of the solar space heating system 
employed in this study appears in Figure 1. The residence 
is heated by forced-hot-air. The flat plate collector 
transforms the absorbed portion of the incident solar 
radiation into thermal energy which is then extracted by the 
working fluid passing over the heated absorbing plate 
(see Figure 2) . This heat transfer is accomplished with 
air, the ultimate delivery medium. When satisfying immediate 
heating requirements of the dwelling, air circulates 
directly between the collector and dwelling - cycle A - if 
useful heat is available. When no solar energy is directly 
available, the room air is circulated through the storage 
medium - cycle B - if useful heat is available there. Should 
there be a requirement for space heating beyond the capacity 
of both these processes, auxiliary heat - cycle C - is 
provided to make up the deficit. Finally, when solar 
energy is available for which there is no immediate need, 
the working fluid is circulated through the storage medium 
transferring thermal energy thereto. As structured, the 
absence of independent heat exchangers (the pebble bed is 
at once the storage medium and a single heat exchanger) 



oua 

mk 
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> residence 



CYCLE 

A heating dwelling with collector 

B heating dwelling with storage (no 

useful heat available from collector 
directly) 

C heating dwelling with auxiliary source 

(not enough useful heat available from 
collector or storage) 

D Heating storage with the collector 

(dwelling requires no heat) 



Schematic diagram of the solar space 
heating system used in this study. 
Operational modes are listed in order 
of precedence. 



FIGURE 1: 
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disallows the ability to charge and discharge storage 
simultaneously as is possible in liquid-type systems. Since 
the transport fluid can also be used as the delivery medium, 
however, no advantage is lost. 

Auxiliary heat, necessary for those periods when no 
direct or stored solar energy is available from the system, 
is provided by a conventional heater. Mounted in the outlet 
duct as shown in Figure 1, this heat exchanger (electrical 
resistance or gas/oil fired) has an auxiliary inlet in order 
that circulation through storage does not occur. An 
exception to this logic could be the case where off-peak 
electrical rates make charging storage with auxiliary, 
off-peak energy economical thus supplementing storage input 
from the collector. 

The storage bed is considered to be located within 
the heated space, therefore convective heat losses from 
the unit are regarded as an uncontrolled gain to the room. 

As shown, the location of the circulation fan would cause 
a reduced-pressure in the collector. Collectors for such 
systems are difficult to seal and this configuration would 
allow infiltration of cold air. However, this has the 
effect of lowering the overall collector temperature 
resulting in more efficient operation at the expense of 
lowering the collector outlet temperature. If one intends 
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to collect the thermal energy at the lowest useful 
temperature anyway, this method would appear to offer an 
advantage. This simulation does not involve such 
considerations, nor does it account for pressure-drops 
through storage or elsewhere. In this regard, it is 
assumed peripheral calculations have addressed compatability 
between mass flow rates used and various flow areas in the 
system (collector and ducting cross-sections, pebble-bed 
effective flow cross-section, etc.). 
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III. SIMULATION OF THE SYSTEM COMPONENTS 



A. Collector 

Several flat-plate collectors are available for air-type 
solar heating systems. This simulation is written to 
accept such units when the collector characteristics have 
been reduced to expressions for U L and F' . The derivation 
of these terms and other factors for a given collector 
geometry is given in [1]. The particular collector 
geometry chosen for this simulation is shown schematically 
in Figure 2. 

To reduce computing costs, the empirical relationship 
developed by Klein for the top loss coefficient, (45°) 
was used: 






N 



1 



-1 



U t (45°) 



(344/T p ) [(T p -T amb )/(N+f)] 



0.31 



*w 



[e +0.0425N(l-e ) ] * + [ (2N+f-l)/e ]-N 

P ir ^ 




( 1 ) 



where 



f = (1.0 - 0.04 h w + 5.0 (10) 4 h w 2 ) (1 + 0.058 N) 



( 2 ) 
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( 2 ) 






e 2 ' h 2' T 2 



U 



b 



U L = U t + U b 



F' = 



U, 



1 + 



h, + 



1 " r l/h 9 + 1/h, 



FIGURE 2: Schematic of Collector Used in This 

Simulation (Elevation) 
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and 



h = 5.7 + 3.8 V 
w 



( 3 ) 



(See NOMENCLATURE section for definitions) 

As this expression corresponds to a collector tilt 
angle, 6, of 45°, an additional equation to interpolate 
for any tilt angle (and, incidentally, for plate emissivities 
other than 0.95) was employed [1]: 

U t (B) = U t (45°) II - (6-45) (0.00259 - 0.00144) e )] (4) 

where f is dimensionless, V is in m/s, all temperatures 

are in °K and 1\ T , U. (45°) and U. (6) are in units of 
W/m 2 -°C. V t • t 

In equation (1) the mean plate temperature, T p , is not 
known. As a result, an iterative solution is used in the 
determination of (45°), (see Appendix). Wind velocity, V, 
is taken as a constant year round but is easily varied if 
desired. The expression used for the back loss coefficient 



Thus the overall collector loss coefficient, U L , is 
determined to be: 



is : 




(5) 




( 6 ) 
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The term h r , representing the radiative heat transfer 
coefficient between the back and absorbing plates (the 
sides of the duct being neglected in the heat transfer) , is 
determined from the righthand-most expression below: 

a(T 1 +T 2 ) (T 1 2 +T 2 2 ) _ 4aT 3 

h = — * (7) 

l/e 1 +l/e 2 -l l/e 1 +l/e 2 -l 

where T, the mean temperature for radiative heat transfer, is 
taken as 340 °K. This assumption is not critical for the 
varying conditions encountered in normal operation [1]. 

The convective heat transfer terms - h-^ and h 2 -, taken 
to be equal, are based on a fully developed forced turbulent air flow 
between the parallel, flat plates. For the mass flow rate and all 
the collector sizes used in this simulation, the Reynold's 
Humber (based on the hydraulic diameter) was greater than 
2100 and the length-to-diameter ratio is larger than 10, 
thereby justifying the assumption. If this is not the case, 
one must consider the effect of the developing boundary 
layer (see Section 4.13 of [1]). The resultant expression 
is : 



h l = h 2 * Nu (* air /D H ) , W/m 2 -°C 



( 8 ) 



-25- 



where 



N u = 0.0158 



m D 



H 



A xs » 



0.8 



(9) 



With these relations and that for U-^, F' is determined: 



F' = 



U. 



1 + 



h l + 1/h- + 1/h 



(1C) 



2 "r 

Depending upon the current temperature of the room or 
storage "cool side" (either of which might be the collector 
inlet temperature depending on the operating mode in effect) 
the energy absorbed by the collector may or may not be 
available at a useful temperature. The expression for useful 
energy gain from the collector per unit area in either case 
is [1] : 



- f r I s - V T c,i - W1 



( 11 ) 



where F R is the collector heat removal factor equal to II] : 



G C 



F R " 



^ 11 - exp (- U t F'/G C )] 



( 12 ) 



U T 



In the simulation, if the derived quantity Q u is negative 
or zero, the collector is not operated. Assuming S, the 
absorbed solar energy (see equ. (24)) is positive the inter- 
pretation of a negative value of Q u is that if run, the collector 
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would be operating at an average temperature too high 
relative to the energy being absorbed to result in an 
increase in fluid temperature. Expressed differently, more 
thermal energy would be lost to the surroundings than is 
being absorbed by the collector. 

Due to the long-term nature of the energy balance 
simulated, the thermal capacitance of the collector was 
neglected. This is justified because the thermal energy 
"lost" in bringing the collector to operating temperature 
in the morning, is subsequently returned to the system at 
the end of the solar collection period. If the time 
response of the system were of interest, for control purposes 
for example, this simplification would not be valid. 

Neglected also was a consideration of the trade-off 
involved between collecting marginal amounts of useful 
energy and the electrical energy expended by the circulation 
fan to do so. 

Collector tilt was taken to be about 20° greater than 
the latitude which emphasizes that this simulation intends 
to size the collector for midwinter heating. Between the 
practical range of tilt angles (latitude + 23.45°), 
particular angles favor particular collection schemes. 

For example, 5 = (latitude + 23.45°) results in the largest 
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ratio of absorbed-to- incident solar energy occuring on 
December 21 - this favors space heating. If $ = (latitude - 
23.45°) , the largest ratio occurs on June 21 and favors 
energy collection on that day. 

This simulation assumes an absorber plate emissivity 
for long wave radiation of 0.95. With selective surfaces 
available at higher cost, this value can be depressed while 
maintaining reasonable solar absorptance values. When 
mass-production techniques are employed to reduce the cost, 
the significant increase in absorbed energy should justify 
the additional expense. 

The effect of dust on the coverplates is accounted 
for as is shading due to the collector frame work. The 
number of covers is taken to be two. Ground reflectance [1] 
is taken as 0.2 year-round which is conservative because 
a value of 0.7 is suggested when there is snow cover. 

See Section III. E . for the equations used to determine 
S, the absorbed portion of the radiant energy SI incident 
upon the outermost coverplate of the titled collector. 
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B. Storage 

The pebble bed storage model is taken from Mumma and 
Marvin [2J. Two difference equations comprise the mathematical 
model with one derived from an air-side energy balance 
and the other from the pebble-side. The equations used 
are [2] : 



T f 

m+1 



(Tf -T fa )exp[-(h v AL)/m C f 
m m 



n) ] + T. 



m 



(13) 



and 



T, (t+At) 
D m 



[A CftT^iTl-Tf^tTD-O^ApCr^d)-^)] 
(pAL/n) (l-e)C fe 



+ T fa (t) 



(14) 



where subscript f refers to the air, subscript b refers to bed 
material (rock) , t is time and subscript m refers to one of 
the n segments into which the storage bed is divided for using 
the model. Other terms are defined in NOMENCLATURE. 



The term h v is the volumetric heat transfer coefficient 
defined by [ 1] : 



h v * 



G/A 



0.7 



W 

m^-°C 



650 



(15) 
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The procedure to use these equations is [2]: 

(a) Assume an initial bed temperature distribution 
at time equal zero. 

(b) For a specified bed inlet fluid temperature, 
evaluate the fluid temperature entering and 
leaving all n sections of the bed using 
equation (13). 

(c) Evaluate the new bed temperatures at each 
section after a time increment At using 
equation (14). 

(d) Step forward in time by Ax. 

(e) Return to step (b) and repeat. 

The "specified bed inlet fluid temperature" is either 
the room temperature or collector outlet temperature 
depending upon the operating mode. The time step. Ax in 
seconds, can be any value between 0 and 3600 in this 
simulation, corresponding to the various operational 
profiles which are possible. In order to insure that the 
solution be stable, the time step in use must satisfy 
the following, [2]: 

(pAL/n) (l-e)c. 

Ax < (IS) 

m c f [1-exp (-h v AL/m d f n)] 



DUG 















* 
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The computer simulation used automatically determines 
the largest number of storage sections, n, which will 
insure a stable solution for the largest possible time 
step which may occur - 3600 sec. 

In a well-designed storage bed the pebble size is small 
enough that the temperature gradient in the individual pebbles 
is insignificant [1], Physically this means the only 
significant resistance to heat transfer to or from the 
storage medium is that due to the convective heat transfer 
phenomenon occuring at the surface of the pebbles. 

The convective heat transfer coefficient associated 
with the surface of the pebbles represents a thermal 
resistance over which we have little control. On the other 
hand, we can manipulate to some extent the thermal response 
of the interior of the pebbles by advantageously choosing 
the size and shape of the pebble in relation to its material 
composition. It is generally accepted [3] that if the Biot 
number (representing the ratio of internal resistance to 
external heat transfer resistance) is less than 0.10/ then 
computational error resulting from assuming a uniform 
temperature within the pebble will be less than 5%. This 
assumption has been made here in order to represent the 
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volume trie heat transfer coefficient, h y , by equation (15) . 

If the pebble size is to be varied to examine the effect, 
this should be taken into account. 

The void ratio, e, is chosen to be 0.3. This corres- 
ponds to a body-centered-cubic packing arrangement of 
exactly similar spheres. Practically, this is the smallest 
ratio which might occur. 

Early trial simulations included unsteady conduction in 
order to destratify the temperatures in the storage segments 
while storage was idle. The standard explicit numerical 
technique used 13] introduced unacceptable errors due to 
the relatively large longitudinal dimension of the storage 
segments. This length is dictated by other stability require- 
ments. An implicit method may have improved the results but, 
no such capability is presently included in this simulation. 
Though less elegant, it is believed no misrepresentation of 
the storage response occurs due to neglecting this effect. 

During the heating season, taken here to be September 
through May inclusive, the storage bed is assumed to be 
located within the residence such that convective losses 
from storage are not "lost” at all but contribute to the 
thermal energy requirements of the dwelling. 
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For the months of June, July and August it is assumed 
these losses take place to the ambient temperature. 

In actual practice the storage bed would probably be 
located in the basement of the dwelling where some of the 
convective loss would be recovered but a significant 
portion would be lost to the ground. In order to avoid 
unnecessary complication, the rather artificial device of 
locating the storage bed based on the season was used in 
order not to unreasonably misrepresent one operating 
mode more than the other. 

All the fluid and storage material (rock) properties 
were considered constant at mean system operating 
temperature values (see Table 2) . 
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C . Residence 

The residence is modeled as a degree-day heat sink. 

It was chosen to correspond in that characteristic, and size, 
to that used by Butz 14] in his water-type system (17,000Btu/ 
°F-day = 374 w-hr/°C-hr) . Assumed to be a box-like structure, 
it has dimensions of 45' x 40' x 8' with standard 
construction materials, 3 inches of insulation, 2 x 4's 
on 16-inch centers and 15% fenestration. The effects of 
heat loss through the floor, infiltration of outside air, 
solar radiation through unshaded windows and internal heat 
generation from occupants and electrical appliances were 
considered by Butz. Considering only the walls, roof and 
windows as media through which conductive and convective 
heat transfer may take place, one arrives at a figure of 
18,200.0 Btu/ °F-day; a 7% error which overpredicts the energy 
need. It appears a much less strenuous calculation may be 
adequate to characterize a residence. 

The model produces data based on the following 
assumptions regarding a desirable home heating scenario: 

(i) During the heating season the comfort range is 
between 67. 0-70.0 *F (19.4 - 21.1 °C) between the hours 
of 0630 and 2230 and 62.0 - 65.0°F (16.7 - 18.3 °C) for 
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the remainder of the 24 hour period. For the purpose of 
determining storage loss to the surroundings or collector 
inlet temperature during direct residence heating from the 
collector, the mean of the desired temperature range in 
effect is used. It is tacitly assumed that the residence 
is within the desired temperature range at all times by 
some combination of direct, stored or auxiliary energy. 

(ii) During June, July and August the desired room 
temperature is between 70.0 - 73.0°F (21.1 - 22.8 °C) 
between hours of 0630 and 2230 and 65.0 - 73.0°F (18.3 - 
22.8°C) for the remainder of the 24 hour period. Although 
during these months storage temperatures are allowed to 
be 180° F (a temperature too high to be acceptable for space 
heating), space heating would in fact be possible with a 
simple mixing damper in the ductwork. The effect would be 
to divert only part of the fixed air flow through storage. 
The present model does not account for this possibility and 
simply provides auxiliary heat when necessary during these 
summer months. 

(iji) If the ambient temperature is greater than the 
desired room temperature but less than the maximum acceptable 
room temperature (73.0°F) no input of any kind is provided. 
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D. Auxiliary Heat 

Auxiliary heat is supplied in the amount necessary 
to make up the deficit between that which the solar heating 
system is able to provide (including uncontrolled storage 
losses to the room) and that required by the residence. 
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E. Insolation and Weather 

The method of Liu and Jordan [5] was used in this 
simulation to provide insolation input. The development 
and application of this method follows directly that of 
Felske [6]. In this approach, actual hourly values of 
insolation for an entire year are unnecessary. Instead, 
this method requires only that the monthly average values 
of H ; and K d be determined for a given location. These 

terms are defined in NOMENCLATURE and [5] explains how 
they are derived. The location used in this study was 
arbitrarily chosen to be New York City (40.77°^). 

With the above parameters determined for each month, 
one is able to calculate for an average day of that month, 
the hourly flux of solar radiation incident upon the 
collector from the equation [6j: 



SI = r d H{ [R b (l-K d /K t ) ] + [1/2 (l+cos3)K d /K t ] 

v ' v ^ / 

beam diffuse sky 

+ [1/2 (l-cose)p ]} 

v 

diffuse ground reflection 



( 17 ) 



where it is assumed that r fc = 



r d and : 
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cos 0 

R b = (18 ) 

sin a 

cos0 t = cos6 cosu) (cosy sin<f> sing + cosiji cosg) 

+ siny sing cos6 sinoj 

+ sin6 ( sin<f> cosg - cosy cost)) sing) (19) 

sina = sin6 sin<J) + cos6 coscf) cosoo (20) 



5 = 23.45 sin [360 (284 + N)/365] 



( 21 ) 



r 



d 



7T COSU) “ COSU) 
S 

24 sinu)„ - u)„cosu) 
s s s 



( 22 ) 



cosu) s = - tan<f> tan6 



(23) 



where N is here the number of consecutive days to mid-month, 
measured from January 1. When the quantity SI is negative 
or zero the collector was not operated. 



Of course the insolation striking the tilted, outermost 
surface of the collector is not the energy which ultimately 
reaches the absorbing plate. The hourly rate of radiation 



-38- 



striking and being absorbed by the collector plate (accounting 
for dust, shading, reflections and absorption) was computed 
from the relation [6]: 



S " (Ta) eff,b r d HtR b (1 -V K t>] 



(24) 



r d H 



+ ( Ta ) e ff,d [ (l+cos8)K d /K t + (l-cos6)p gr ] (1-s) (1- 



d) 



where from (6]: 



(ta) 



eff 



T(6 t ) a(0 t ) 

l-[l-a(0 t ) ] p (60° ) 



(25) 



and 



Td t ) = 



1-P(0 t ) 



l+(2N-l)p(0 t ) 



(26) 



1 


sin 2 


<J> 

to 

1 

CD 

ft 


tan 2 


<e 2 -e t > 


2 


sin 2 


<w 


tan 2 


(e 2 + e t ) 



(27) 



©2 = sin ^ 



sin0. 



N, 



(28) 



in which N equals the number of glass cover plates and other 
terms are defined in NOMENCLATURE. 
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Reference [6] indicates that the directional absorptivity 
of the collector surface, a(.0^.), can be represented by a 
curve for a typical surface. The equation provided by Felske 
to represent this typical curve is: 

20.03 

a(0.) = 1.0 t— (29) 

(100.03 - 0 t )- L * J 

where 0^ assumes its actual hourly value for beam radiation 
but is assigned the value of 60° when used to represent the 
diffuse radiation component [6]. 

The hourly ambient temperature distribution is more 
random than insolation and it can probably be argued that 
one reasonable distribution of the average daily 
temperature is about as good as another. In this simulation, 
the scheme illustrated in Figure 3 was used. The average 
daytime temperature and nighttime temperature were cast 
on a sine wave peaking at 1300 hours (Note: throughout this 

simulation, hourly values are referenced to solar time) . 
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tewperature 




FIGURE 3: 



Distribution of Ambient • 
Temperature Through the 
Day 
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TABLE 2 

Simulation Design Parameters 



c b 


0.837 KJ/Kg-°C 


c f 


1.012 KJ/Kg- °C 


d 


0.02 


Insulation thickness, collector 


0.076m 


k air 


0.029 W/m- ° C 


k ins 


0.043 W/m-°C 


m 


1.2 Kg/s 


D 


0 . 0508m 


N 


2 

(number of covers) 


n g 

Plate spacing, collector 


1.526 

0.01m 


s 


0.03 


V 


5 . 0 m/s 


a (60 0 ) 


0.8346 


3 


60° 


e 


0.3 


e g 


0.88 


e l 


0.95 


Y 


0.0 
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p (absolute viscosity) 



1.912 (10)" 5 Kg-m/s 



* 

P 



40 . 77° N 
2400 Kg/m 3 



P 



gr 



0.2 
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TABLE 3 



Approximate Time Constants for Major System Components 

e 

Values for collector and storage are based on the basic 
component parameters, assumed constant at approximately the 
mean operating temperature associated with each component. 

The value for the residence is calculated based on 
the time to change the interior temperature of the author's 
home a given number of degrees over a period of time with 
constant ambient temperature and no heat input. 



Collector: 



1/4 hour 



Storage Bed 



1 hour 



Residence: 



10 hours 



- 44 - 



IV. OPERATIONAL LOGIC AND STRATEGY 

The flow chart of Figure 4 represents the operation of the 
solar system. A breakdown of the computer logic used to 
implement this basic logic appears in the Appendix. 

This simulation operates by accounting for the following 
heat transfers from and into the residence unit: 

a. The heat loss from dwelling driven by the 
difference between the desired room temp- 
erature and ambient temperature. 

b. Heat supplied directly to the dwelling from 
the collector on an as-needed basis when 
available. 

c. Heat supplied to the dwelling from the storage 
bed in the form of either uncontrolled con- 
vective losses or forced-hot-air plus the 
uncontrolled losses. 

d. Heat supplied from the collector to the 
storage bed. 

e. Auxiliary heat as necessary to meet what 
remains of the heat load should sources 
(b) and (c) prove insufficient to meet the 



entire demand. 
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figure 4 Flow Diagram of Overall Simulation Logic 
(Refer also to Figure A. 2) 
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f. Total useful thermal energy obtaining from 
the collector for all purposes. 

The only temperatures monitored or used as controlling 
input are those of the storage bed segments. The temperature 
of the residence is assumed at all times to be within the 
comfort range in effect for a given season. 

To implement the insolation representation of Liu and 
Jordan I 5 ] , it has been assumed that the monthly average 
of the hourly insolation and ambient temperature applied for 
one 24-hour period predicts the system behavior for the 
entire month when the results of this "average” day are 
multiplied by the number of days in the month. Butz [ 4 ] 
used a similar extrapolation with hourly Weather Bureau data. 

A discontinuity occurs in the system operation in going 
from one month into another due to the change in these monthly 
average values. Furthermore, the final temperature profile 
in the storage bed carries into the succeeding month. For 
these reasons the first days of operation in any given month 
do not represent a true average daily response. If the 
simulation is operated for several days utilizing the 
insolation and ambient temperature profiles of that month. 



the influence of the starting conditions is eliminated. A 
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measure of the degree to which the starting conditions have 
been obscured is the correspondence between storage bed 
temperature profiles separated by 24 hours of operation. 

For the ranges of parameters and values of constants 
used here, 4 to 5 days operation served to damp out starting 
effects in most cases. Therefore, the following strategy 
was used: 

i. For each month, 9 full days of hourly 
operation were simulated. 

ii. The last run (corresponding to the 9th 
day) is used as the "average" day. The 
results obtained for this day are multiplied 
by the number of days remaining in the month. 

iii. The results for the entire month are obtained 
by adding the cumulative results of the first 
9 days to those obtained in (ii) . 

During the heating season, the maximum allowable average 
temperature of the storage bed is chosen to be 145°F. In 
the limit, then, this is the highest possible forced-air 
temperature which can be derived from storage. This temp- 
erature is at least 20°F below that which outlets from a 
conventional forced-hot-air heating unit where common 
practice indicates a rise in temperature of 100 °F or more 
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through the furnace 18 J . This lower operating temperature 
allows more efficient operation of the collector but requires 
a higher-than-conventional mass flow rate of warm air into 
the room. The codes 18] restrict the rate of room air change- 
out and must be consulted. If upon supplying the hourly 
room needs the collector is able to bring storage to this 
temperature, the collector ceases to operate. Subsequent 
collector-to-storage operation is only to the extent necessary 
to maintain this temperature in the face of demands made upon 
storage. This mode of operation causes the low usefulness 
ratios (useful heat collected/incident thermal energy) in 
early fall and late spring. 

During the summer months when this simulation does not 
supply solar space heating, the maximum storage temperature 
is allowed to be 180 "F for the purpose of domestic water 
heating. Higher temperatures are possible if pressure relief 
and mixing values are used in the water system but then the 
collector would operate less efficiently. 180 °F was chosen 

to avoid unnecessary complication. No attempt is made to 
model a realistic daily hot water usage scenario in the summer. 
Instead, the heating load presented to the bed is that of an 
increased convective heat loss coefficient corresponding 
approximately to an 80 gal/day draw-down of 100 °F water spread 
evenly through the 24-hour day, added to the previously used 
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coefficient. Accurate assessment of the water heating mode 
does not contribute to the objective of this study; this 
portion of the simulation is included only for continuity. 

If the collector can heat the dwelling directly but is 
unable to meet the full requirement, auxiliary heat is used 
to make up the deficit. In this situation the only contri- 
bution from storage is the uncontrolled convective loss. 
This strategy is used in order to obtain all the useful 
solar energy available and does not consider the electrical 
energy required by the circulation fan to obtain this 
marginal amount of thermal energy. 

Several arbitrary choices have been made as to what 
constitutes a desirable system operation. 

i. If when heating storage with the collector the 
increase in average bed temperature corresponds 
to a thermal energy input of less than 3.3 W-hr 
in one minute (an hourly rate of 200 W) , the 
system is not operated. 

ii. When heating the room from storage a 1% overshoot 
of the room's needs is tolerated. The overshoot 
shows up as excess heat which is subtracted from 
the otherwise useful energy gain of the system 
and results from the time step used (see iv) . 
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iii. If at any time while heating the room from 
storage the thermal energy added in one minute 
is less than 5% of l/60th the total hourly room 
need, then the system is relaxed and the 
auxiliary source provides the remainder of 

the heat load. 

iv. All heat transfers to and from storage are 
accomplished in 1-minute intervals. If some 
fraction of a minute over the hour is used or 
an unacceptably large overshoot in energy input 
occurs, the 1-minute interval in which this 
occurred is nullified and the process continues 
in 5-second intervals. Any further overshoot 
resulting during such a reduced interval is 
accepted and appears as an accountable excess. 

When heating storage from the collector, the inlet 
temperature to the collector is held fixed for the duration 
of the hour at its value at the beginning of the hour. 

In fact this temperature increases as the bed is heated 
and, if strictly accounted for, could result in a decrease 
in useable thermal energy available to storage. Errors due 



to this assumption are considered small considering the 
small change in temperature of the storage bed which occurs 
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in one hour. So, in the interest of decreasing computation 
time, the effect was neglected. 



V. RESULTS 



Figure 5 demonstrates the ability of various collector 
sizes to meet the heating demands of the residence used 
in this simulation. As collector size increases a greater 
proportion of the room needs can be accommodated. If 
analysis of cost tradeoffs and regional weather character- 
istics indicate that it is advantageous to use solar power 

for approximately half the residential heating needs in 

. 2 
mid-winter, then this figure indicates that the 60 m 

collector should be employed. 

The curves of the ratio of solar energy obtained to 
room needs for various collector sizes were derived using 
a fixed storage volume. Figure 6 indicates, however, that 
mid-winter system capability for various collectors is 
relatively insensitive to the size of storage. The reason 
for this result is that for the collector sizes examined, 
no collector provides on a daily basis more energy than 
the residence consumes in that time. There is then no 
net accumulation of energy in storage and the total of the 
daily solar energy obtained reflects the room needs only. 

The uppermost, flat curve of Figure 5 represents 
the practical limit of system output. Achieving this 
output requires a collector of such size that storage is 
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FIGURE 6 Daily Integrated Energy Ratios, Averaged by Months for Fixed Collector Area 
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brought to a fully charged state each day while the 
immediate daytime dwelling heat load is met as well. A 
"fully charged" storage bed is one at the maximum allowable 
average temperature. It is assumed in presenting this 
limiting case that sufficient storage volume has been 
assigned to the system to meet, in full, the overnight 
heat load if the bed is initially fully charged. 

In Figure 5 the QU/QROOM curve (see SIMULATION 
NOMENCLATURE for definition of terms) for the largest 
collector used is seen to behave differently than the 
other two curves during April and May. The "cross over" 
observed results when the storage bed approaches the 
maximum average temperature allowed. When this occurs, 
very little further solar energy input to the system takes 
place. 

A QU/QROOM ratio exceeding the maximum practical 
value of 1 is seen to occur in late spring for all three 
collector sizes examined. This is because in those months 
enough solar energy can be obtained from the collector to 
meet the room needs and allow a net accumulation of energy 
in storage. Note that as collector size increases, the 
QU/QROOM curve approaches the limiting curve. 
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The curves of QU/QI of Figures 5 and 6 represent the 
ratio of solar energy collected to that incident upon the 
tilted collector surface. This ratio is a measure of collec- 
tion efficiency. For fixed system mass flow rate of air and 
storage size, increasing collector area decreases collection 
efficiency because of the higher mean temperature of collec- 
tion and storage (even if panels are added in parallel) . 

Efficiency is generally lower in early fall and late 
spring because the higher insolation and storage temperature 
approaching saturation cause a higher mean temperature of 
collection. This, in turn, drives a higher loss rate to the 
surroundings despite relatively warm ambient temperatures. 

In late fall and early spring the relation between ambient 
temperature, available insolation and the lower average temp- 
erature in storage combine advantageously to produce the 
highest efficiency. In mid-winter the low ambient temper- 
ature drives large collector losses despite the rather low 
mean temperature of collection (resulting from the lowest 
insolation and storage temperature levels of the year) . 

In Figure 6, the higher efficiency of larger storage beds 
in spring reflects the increasing level of collector-to-storage 
heating activity and lower collection temperatures which result 
from distributing about the same amount of collected energy 
throughout a larger mass of rock. The approach saturation which 
will occur as storage size decreases, rapidly increases the 
efficiency disparity. 
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Figures 7 and 8 illustrate the monthly totals of 
hourly solar energy supplied by the system in excess of 
that required. During early spring and late fall the 
excess is primarily caused by the uncontrolled convective 
storage loss. During the winter when all the storage loss 
is of use to the dwelling, any excess occuring is caused 
by an "overshoot" of storage energy input to the room 
during active heating with storage. This overshoot results 
from the time increment used in the simulation (see pg. 50 
Sec. IV) and is not a characteristic of solar systems. 

The excess in spring and fall is a real effect if 
storage is contained within the dwelling. These losses 
are annoying and are counterproductive to any air- 
conditioning employed during the day. 

All three curves of Figures .7 and .8 extend to about 

5 

1.4(10) W-HR in September, off the scale presented. This 
occurs because storage was taken to be fully charged at 
the commencement of this month. The high convective heat 
losses which are a consequence cause the high excess 
energy value. 

Such excesses do not occur in May because storage 
temperature levels do not naturally reach saturation 
with the available solar energy and room requirements 



Collector Dimensions: 3m x 10m 
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FIGURE 7 Monthly Quantity of Excess Energy, Average by Months, for Fixed Storage Volume 



Storage Dimensions: 2.0m x 6,0m 2 
Storage Dimensions: 2.5m x 6-8m^ 
Storage Dimensions: 3,2m x 7,5m 2 
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FIGURE 8 Monthly Quantity of Excess Energy, Average by Months, for 5ixed Collector Area 
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are higher in May than September. Starting the month of 
September with a fully charged storage bed simulates having 
used storage during the summer for domestic water heating 
at the maximum allowable average bed temperature of 180 °F 

The two larger collectors are seen to provide similar 
amounts of energy in excess of that needed for the month of 
May. This results because, again, the storage bed is 
approaching the saturation level. Once near or at this 
level, a larger collector will not result in more storage 
losses. 

In Figure 8, an apparent anomoly occurs in late spring 
when the excess energy for the largest storage unit falls 
below that of the other two units. An examination of the 
hour-by-hour energy transfers to the room from storage 
indicates that the "overshoot" flaw is the cause (see pg. 50, 
Sec. iv) . A lesser percentage of excess energy input to 
the dwelling occurs each hour with the largest bed compared 
to the two smaller units. This behavior is evidence of the 
lower overall storage temperature as bed size increases, all 
else being equal. 

The arbitrary choice of storage size used in this 
simulation happens to span the range over which an inversion 
occurs in excess energy output from storage. It is observed 
in Figure 8 that from October through December the excess 
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energy first increases with increasing storage volume then 
decreases. Again, inspection of hour-by-hour energy 
transfers indicates the cause. The smallest bed can provide, 
in the first several hours after sundown, more hourly energy 
output than larger units which have a lower average temper- 
ature. This is so despite the fact that all beds contain 
about the same total stored energy. 

The higher hourly output from the smallest storage bed 
results in more hourly overshoot. But since this bed is 
depleted sooner, it is less often able to provide heat. The 
cumulative overshoot is consequently less than for larger 
units. 

The largest unit operates to heat the room the most 
number of times, but each time the overshoot is small because 
of the relatively lower temperature of the energy provided. 
Cumulatively, less total overshoot occurs than with the 
mid-size bed but more than that occurring with the smallest. 

Figure 9 shows the degree in which room energy needs are 
met with the desired collector and various storage volumes. 
This figure applies to the New York area for September and 
October starting with a storage bed at room temperature. 

For the collector size and residence of this example, the 
only months during which more solar energy is provided than 



Fixed Collector Area: 60m 
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FIGURE 9 Ability of Various Storage-Collector Combinations to Meet Room Heating Needs 
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that required to meet the daily room needs is September and 
May. Either of these months may be used to "sense" the 
presence of the storage bed, no longer masked by the over- 
powering effect of room needs. Relative stroage size plays 
very little role during mid-winter because then storage only 
serves to distribute solar energy stored during the day 
through part of the immediately following night with no 
carry-over of energy into the following day. It is clear 
that some storage volume is necessary but that size has yet 
to be determined. 

The dashed, vertical line of Figure 9 indicates the 
collector-to-storage ratio for which no more solar energy is 
collected daily by the system than is required for an 
average day in September. Allowing nothing for "cloudy 
days" carry-over storage capability, this system should have 
a storage volume in September equal to the collector area 
times the factor 1/35 (the dimensions being in meters) . 

A storage sizing criterion which optimizes mid-winter 
use of the solar heating system would seem to be a more 
rational criterion than the above. It has been shown, 
however, that mid-winter system operation is relatively 
insensitive to storage size. Within the range of storage 
sizes examined, the ratio of storage output to room needs 
increases less than 1% (see Table 4) when storage volume increases 
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twofold. So, such a criterion is elusive at best. Intuiti- 
tively, there are limits to the range of useful storage 
temperatures and this suggests an upper and lower bound to 
storage volume exists. 

If an upper temperature limit is placed on the bed 
(as dictated by the codes as regards delivery temperature 
or room air change-out rate, etc. : in this simulation 

145” F was used) there is a lower storage volume which 
will attain, on average, this saturation temperature using 
all the energy available to storage on a typical mid-winter 
day. However, because the collector operates inefficiently 
due to the high temperature it receives from storage, the 
total energy stored is less than could have been obtained 
for any larger storage bed given the same conditions. 

The lower limit on storage bed temperature is that 
which provides a comfortable influx of heating air and 
depends on relative humidity. Generally, normal body 
temperature represents a good estimate of this lower limit 
on storage temperature. The storage volume for which this 
temperature obtains when all the energy available to storage 
on a mid-winter day is accepted, is the upper limit on 
storage volume. This provides the most efficient system 
operation based on collected energy but the energy stored 
at this temperature is only marginally useful. 
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Lacking a definitive raid-winter storage-sizing method, 
it is suggested that the "September storage sizing criterion" 
above be adopted as a method of determining storage volume 
given a collector area. 

In this simulation, using a 5m x 12m collector the 

. 3 

storage volume thus indicated is 1.7m . If it is assumed 

that the energy available to storage in January for a 
2 

6.6m x 2.5m storage volume (see Table A.l) remains 

3 

constant as the storage volume decreases to 1.7m , the 

o 

resulting average storage bed temperature attains 115 F 
(if it starts at room temperature) . Of course, due to 
collector inefficiency the collected energy would be less 
as would be the temperature attained. 

Thus, an even smaller storage bed could be used with- 
out saturation occurring but less total energy would be 



stored. 
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TABLE 4 



Energy input to room from storage, between sunrise 
and sunset, as a percentage of room needs during 
that period. Values are for January with 
combinations of storage and collector sizes as 
indicated. 



3 

Storage Volume (m ) 
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00 
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VI. CONCLUSIONS 

From Figures 5 and 6 it can be concluded that collector 
area requirements to meet a known energy need can be deter- 
mined without regard to storage size. 

It is unnecessary to consider increasing storage 
size to allow for the "carry over" of solar energy into 
subsequent "cloudy days" because for economically sized 
collector areas designed for mid-winter needs there is, 
on the average, no energy available for subsequent days. 

Only if more than 100% of the daily room requirements 
can be supplied from the system in mid-winter is there, 
on the average, a "carry over" capability which can be 
enhanced by increasing storage size. Collector size is 
the single most important consideration in installing 
this solar heating system. 

Once the collector area necessary to meet the energy 
need has been determined a compatible storage volume must 
be derived. "Compatible storage volume" is taken to be 
that volume which results in a QU/QROOM ratio of 1 in 
September (see Figure 9) . It must be realized that such 
compatability of collector and storage applies only to 
a system providing approximately 50% of the mid-winter 
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space heating requirement and the location for which this 
analysis was carried out, that is, New York. 

Within the above conditions, it is suggested that the 
ratio 35:1 be used to relate collector area to storage 
volume (dimensions in meters) for a hot air system with 
pebble bed storage. Researchers generally quote suggested 
ratios of collector size per storage volume (or per person 
for strictly domestic water heating applications) when 
speaking of water-type systems. When doing so, the 
validity of the ratio is predicated on specifying the 
location and percentage to be met of the total demand. 

It is for this reason, with a similar claim of validity, 
that the ratio above is presented. 



APPENDIX 



The flow chart included here illustrates the simulation 
of the solar heating system used in this study. 

The numbers appearing in the chart are keyed to the 
various "go to" statements used in the computer program and 
are included in order to guide the reader through the 
computer listing appended. 

The terms used herein are defined in the SIMULATION 
NOMENCLATURE section of the appendix. 

Two pages of output are presented to illustrate their 



form. 
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e = error 




FIGURE A. 1 . Determining the Mean Plate Temperature [lj 
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GRMR4S-QRMFNS+ QPD 

_ . . qmtx^rooh-qrkfhs-qrmfhc 

FIGURE A. 2 (cont.) | qxces-Iqaux) if Q UAX<0 
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TABLE A. 1 

Sample of Yearly Simulation Output for an 
Average Day in Each Month 



Collector Area: 5m x 12m 

2 

Storage Volume: 6.8m x 2.5m 

Temperatures in °C 
Hours are Solar Time 



HOURLY HOURLY HOURLY HOURLY HOURLY HOURLY 

HOURLY HOURLY ROOT. COLLECTED STORAGE COLLECTED AUXILIARY STORAGE 

AMMENT INCIDENT USEFUL HEAT HEAT TO HEAT TO HEAT TO HEAT TO LOSS 

TI"'!' HEAT HEAT MEEDS ROO* SURR STORAGE ROOM TO SURR 

HOUR (C) (H-HR) (W-H») (W-HR) (W-HR) (W-HR) (S-HR) (W-HR) (U-HR) 
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COMPUTER NOMENCLATURE 
FOR FIGURE A. 2 


AREAC 


surface area of collector 


DEGDAY 


heat load of dwelling; (W-HR/°C-HR) 


AT 


difference in temperature between storage bed 
surroundings and average temperature in bed 


At 

s tor 


change in temperature of storage due to 
heating bed with the collector 


aqpb 


heat transferred from storage while storage 
is idle (convective losses) 


AQRB 


heat transeferred to storage from collector 


DLTIME 


time, in seconds, regaining in any given hour 
for performing any of the several modes of 
heat transfer possible in this simulation 


DLTIMH 


one hour in seconds 


IFLAG 


when equal to 1, indicates the collector 
is presently heating storage after having 
heated the dwelling; used as an internal 
logic signal only 


QADINT 


heat transferred in an interval (either 60 sec 
or 5 sec) when storage used to heat dwelling 


QADTOT 


heat transferred while heating room from 
storage 


QAUX 


auxiliary heat required by dwelling 


QAXQRM 


monthly ratio of auxiliary heat supplied to 
room needs 


QI 


solar radiation incident upon collector 


QLOSSH 


Total amount of heat lost from storage to 
room vhiie storage being heated by collector 
(QLOSS refers to the amount lost in either 
60 sec or 5 sec) 



QRATIO 


ratio of energy supplied from storage in 
small interval to the room needs in that 
interval 


QRMFMC 


heat transferred to the dwelling from collector 


QRMFMS 


heat transferred to the dwelling from storage 


QROOM 


hourly room energy requirement 


QROOM1 


hourly room energy requirement accounting 
for the uncontrolled storage loss 


QSLRQR 


monthly ratio of solar energy supplied to 
room needs 


QSTFMC 


heat transferred to storage from collector 


QSUBLS 


uncontrolled convective losses from storage 


QU 


hourly thermal energy extracted from the 
collector by the system = QRMFMC + QSTFMC 


QUQI 


ratio of QU to QI 


QUQRM 


ratio of QU to QROOM 


QXCES 


the hourly amount of energy supplied to the 
residence for which there is no need 
(convective losses or "overshoot" from 
storage) 


QXCSMT 


monthly total of QXCES 


SI 


flux of solar radiation incident upon the 
collector surface area 


SURFAR 


surface area of storage bed through which 
convective loss takes place 


TAMB 


ambient temperature 


TAVCHK 


maximum average storage temperature allowed 


TREF 

TSAVG 


average value of allowable room temperature 
range 

average instantaneous temperature in storage bed 
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